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In  this  paper,  structural  configuration  of  layered  Li-ion  battery  electrode  plates  is  evaluated  by  analytically 
formulating  the  diffusion  induced  stress.  Both  symmetric  electrode  and  asymmetric  bilayer  electrode  are 
discussed.  The  thickness  ratio  and  the  modulus  ratio  of  current  collector  to  active  plate  are  analytically 
identified  to  be  the  important  influence  parameters  on  the  stress.  Applying  a  material  with  smaller 
elastic  modulus  for  current  collector  could  reduce  the  peak  stresses  in  both  current  collector  and  active 
plate.  Increasing  the  thickness  of  current  collector  would  reduce  the  stress  in  itself  while  promote  the 
stress  in  active  plate.  Therefore,  from  mechanical  viewpoint  on  designing  an  electrode,  the  material  for 
current  collector  should  be  as  soft  and  flexible  as  possible.  And  the  thickness  of  current  collector  should 
be  in  an  appropriate  range.  Basically,  it  should  be  as  small  as  possible  on  the  precondition  that  the 
mechanical  strength  is  satisfied.  Finally,  effects  of  three  charging  conditions,  i.e.  uniform,  galvanostatic, 
and  potentiostatic,  on  the  diffusion  induced  stress  is  discussed.  It  is  found  the  maximum  stresses  for 
three  cases  are  linear  to  the  total  amount  of  intercalated  lithium  ions.  Based  on  the  stresses,  an  optimized 
charging  operation,  i.e.  first  galvanostatic  followed  by  potentiostatic,  is  suggested. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Internal  stress  in  Li-ion  battery  electrode  generated  by  interca¬ 
lation  and  deintercalation  of  Li  ions  into  and  from  active  materials 
has  received  considerable  attention  because  it  is  a  direct  cause  of 
cracking  which  leads  to  the  loss  of  capacity  and  the  deterioration 
of  cyclic  performance.  In  the  literature,  the  insertion/extraction 
of  lithium  ions  and  the  induced  stress  in  an  electrode  have  often 
been  modeled  by  using  the  continuum  diffusion  model  coupled 
with  continuum  mechanics  model.  Some  considered  one-way  cou¬ 
pling  models  where  the  diffusion  is  not  affected  by  stress  whereas 
the  others  adopted  more  sophisticated  two-way  coupling  models 
where  stress  and  diffusion  are  mutually  affected.  Huggins  and  Nix 
[1]  presented  a  simple  one-dimensional  beam  model  for  analy¬ 
sis  of  cracking  caused  by  swelling  transformation  strain  in  bilayer 
plate  structures.  Garcia  et  al.  [2]  performed  a  two  dimensional 
numerical  analysis  of  diffusion  of  Li-ions  in  both  electrolyte  and 
spherical  active  materials  in  cathode  with  results  indicating  that 
the  large  stresses  develop  during  fast  discharge.  Christensen  and 
Newman  [3]  developed  a  strongly  coupled  diffusion-stress  model 
to  calculate  stress  in  spherical  active  particles  emphasizing  on  the 
pressure-diffusion  effect  and  the  high  Li-ions  concentration.  Zhang 
et  al.  [4]  conducted  one-dimensional  modeling  of  a  spherical  active 
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particles  and  three-dimensional  simulations  of  ellipsoidal  parti¬ 
cles  under  galvanostatic  Li  ions  intercalation;  the  numerical  results 
for  LiMn204  system  suggest  that  particles  with  smaller  sizes  and 
larger  aspect  ratios  lead  to  smaller  intercalation-induced  stresses. 
Cheng  and  Verbrugge  [5]  have  developed  analytical  expressions  for 
the  evolution  of  stress  and  strain  energy  within  a  spherical  parti¬ 
cle  under  both  galvanostatic  and  potentiostatic  conditions.  Using 
a  diffusion-elasticity  model  with  Butler-Volmer  surface  kinetics 
Golmon  et  al.  [6]  simulated  the  insertion  of  lithium  into  spherical 
silicon  particles  with  results  suggesting  that  the  stresses  depends 
strongly  on  the  reaction  kinetics  and  discharge  rate.  Using  the 
cohesive  model  of  crack  Bhandakkar  and  Gao  [7]  predicted  a 
critical  characteristic  dimension  for  crack  nucleation  in  an  ini¬ 
tially  crack-free  strip  electrode  under  galvanostatic  intercalation 
and  deintercalation  processes.  Qi  and  Harris  [8]  demonstrated 
graphically  the  lithium  concentration  spatial  maps  via  the  color 
variation  in  graphite  and  carried  out  an  in  situ  observation  of 
strain  held  during  lithiation  of  a  graphite  electrode.  Using  the  the¬ 
ory  of  diffusion-induced  stress  and  the  energy  principle,  Yang  [9] 
derived  an  analytical  relation  between  the  critical  concentration  of 
solute  atoms  and  average  damage  size  for  the  insertion-induced 
cracking  and  buckling  in  an  elastic  him.  Seo  et  al.  [10]  performed 
a  hnite  element  analysis  of  intercalation-induced  stress  in  the 
LiMn204  particle  with  the  surface  morphology  reconstructed  using 
atomic  force  microscopy  (AFM),  showing  that  Mises  stress  at  the 
sharply  dented  region  is  an  order-of-magnitude  higher  than  the 
stress  in  smooth  particle.  Zhao  et  al.  [11]  developed  a  coupled 
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diffusion-large  plastic  deformation  theory  with  application  to  anal¬ 
ysis  of  the  induced  stress  and  plastic  deformation  in  a  spherical 
particle. 

The  most  works  mentioned  above  concerned  the  diffusion 
induced  stress  at  particle  level.  Li  diffusion  and  diffusion-induced 
stress  at  the  higher  scale  of  layered  structure  are  usually  neglected. 
It  is  very  important  to  discuss  the  layered  structural  configura¬ 
tion  where  the  multilayer  electrode  plate  includes  active  plates  as 
well  as  current  collector.  Furthermore,  Li  concentration  gradient 
at  the  layer  level  in  active  plate  during  charging  was  experimen¬ 
tally  observed  [12],  indicating  that  all  active  particles  in  the  active 
plate  are  not  charged  simultaneously  in  some  cases.  The  layered 
geometry  is  one  of  promising  structural  configurations  for  newly 
developed  nanoelectrodes  [13].  Therefore,  focusing  on  the  multi¬ 
layer  electrodes  the  purposes  of  this  work  are  to:  (1 )  introduce  the 
spatial  distribution  of  Li  concentration  and  develop  a  general  and 
yet  analytical  solution  for  diffusion-induced  stress  in  multilayer 
electrode,  (2)  identify  the  parameters  of  layered  structure  which 
affect  the  diffusion-induced  stress,  (3)  address  the  role  of  current 
collector  on  the  induced  stress  in  electrode. 


2.  Analysis 


Consider  an  electrode  plate  in  which  two  active  material  plates 
are  bonded  to  a  central  current  collector,  Fig.  1,  where  hi,  h2  and 
hc  denote  the  thickness  of  two  active  plates  and  current  collec¬ 
tor,  respectively.  The  active  plates  are  considered  as  composites 
containing  both  active  particles  and  pore-filling  electrolyte.  Let  the 
thickness  direction  be  aligned  with  the  z-axis  and  the  in-plane  of 
plate  with  x-  and  y-axes.  Each  individual  layer  in  the  electrode  plate 
is  assumed  to  be  an  isotropic  and  elastic  material.  Lithium  ions  are 
allowed  to  be  inserted  into  and  extracted  out  of  the  electrode  plate 
from  one  side  face  of  each  active  material  plate,  Fig.  1. 

The  diffusion  of  lithium  ions  in  the  thickness  direction  of  elec¬ 
trode  is  assumed  to  be  governed  by  Fields  law 


—  -  =  0 
dt  dz 2 


(i) 


where  c  is  the  molar  concentration  of  lithium  ion,  D  is  the  effective 
diffusivity  of  Li  in  the  composite  active  material.  It  describes  the 
macroscopic  transportation  process  including  diffusion  and  migra¬ 
tion  of  lithium  ions  through  the  liquid  phase,  solid  phase  and  the 
solid-liquid  interfaces  in  the  active  plate.  The  initial  solute  concen¬ 
tration  in  the  electrode  is  assumed  to  be  zero,  i.e. 


active  active 

plate  collector 


Fig.  1.  Configuration  of  the  electrode  in  which  two  active  plates  are  bonded  to  the 
current  collector.  Here,  hi  =  h2  in  symmetric  electrode,  while  hi  ^  /i2  in  asymmetric 
case. 

follows  [14] 


c  =  0  for  t  =  0  (2) 

The  boundary  conditions  for  active  Plate  1  in  the  electrode  con¬ 
figuration  shown  in  Fig.  1  are  given  by 

forz  =  h\  (3) 

forz  =  0 

for  the  galvanostatic  charging,  where  in  is  the  surface  current  den¬ 
sity  and  F=  96485.3  C  mol-1  is  Faraday’s  constant,  and 


D—  —  in 
dz  F 


4-o 


c  =  c0  forz  =  hi 


(4) 


D—  =  0  forz  =  0 
dz 

for  potentiostatic  charging,  where  Co  is  the  constant  surface  con¬ 
centration.  The  solution  for  the  ion  concentration  can  be  found  as 


c(z,  t)  = 


inh  i 
FD 


Dt  3 z2  -  h2 

hi2  6 h2 


2 

7T2 


n= 1 


nnz 

x  cos  — —  exp 
hi 


n27T2Dt\  \ 

hi2  J  J 


(5) 


during  galvanostatic  charging,  and 


c(z,  t)  =  c0-  coy 


4(-l)m -1 
(2m  -  1  )tt 


cos 


/(2m  -  1  )ttz\ 

V  2hi  J 


x  exp  - 


(2m-  l)27T2Dt 
4hi2 


(6) 


during  potentiostatic  charging.  The  concentration  in  the  active 
Plate  2  can  be  obtained  from  above  equations  by  replacing  z- 
coordinate  with  (-z-hc). 

The  properties  of  active  material  such  as  Young’s  modulus  could 
be  significantly  affected  by  the  solute  concentration  [  1 5  ].  Therefore, 
for  generality  here  we  assume  that  Young’s  modulus  E,  Poisson’s 
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ratio  v  and  partial  molar  volume  £2  depend  upon  the  concentration, 
in  turn,  are  implicit  functions  of  z. 

In  large  electrode  plate  the  equal  biaxial  stress,  oxx(z)  =  oyy(z) 
which  are  functions  of  z  only,  is  assumed  to  be  generated  by  the 
inhomogeneous  distribution  of  concentration  given  by  Eq.  (5)  or 
(6).  The  other  four  stress  components  are  zero.  It  is  easy  to  ver¬ 
ify  that  the  equations  of  mechanical  equilibrium,  3aJl/9xJ  =  0,  can 
be  satisfied  automatically.  Furthermore,  it  is  suggested  that  from 
the  constitutive  equations  the  two  in-plane  normal  strain  compo¬ 
nents  are  also  equal,  i.e.  £**  =  eyy.  The  six  deformation  compatibility 
equations  lead  to: 


d2Sx> 

dz2 


=  0 


(7) 


This  suggests  that  the  in-plane  normal  strain  varies  linearly  in 
thickness  direction, 


£xx  —  £yy  —  Sq  +  ZK 


(8) 


where  £0  is  considered  as  in-plane  strain  at  the  plane  of  z  =  0  and  k  is 
the  curvature  of  deformed  electrode.  This  suggests  that  the  active 
plates  and  current  collector  plate  have  the  same  bending  curva¬ 
ture  when  they  are  perfectly  bonded.  The  constitutive  equation  for 
biaxial  stress  is 


&xx  —  Gyy  —  E'{Sq  +  KZ )  —  —E'£2c 


(9) 


where  E'  =  E/(  1  -  v)  is  the  biaxial  modulus,  Q  is  the  partial  molar 
volume  of  solute.  The  concentration  related  term  -EQcj 3  rep¬ 
resents  the  diffusion  induced  stress,  written  in  an  analogue  form 
of  thermal  expansion  stress.  This  term  should  be  removed  for  the 
current  collector. 

The  traction-free  boundary  conditions  on  two  side  faces  have 
been  satisfied  automatically.  Since  we  focus  on  the  diffusion 
induced  stress  which  does  not  include  the  effect  of  external 
mechanical  loading  the  end  of  plate  is  mechanically  free  and  the 
mechanical  boundary  conditions  could  be  specified  as 


axxdz  =  0 


OxxZ  dz  =  0 


(10) 

(11) 


Here  we  emphasize  that  the  above  obtained  solutions  for  stress 
and  strain  are  exact  because  all  mechanical  field  equations  of 
three-dimensional  elasticity,  interface  continuity  conditions  and 
boundary  conditions  have  been  satisfied. 


3.  Symmetric  electrode  plates 


In  this  section  we  exam  the  electrode  which  includes  two  active 
material  plates  with  equal  thickness  by  setting  h2  =  h\.  The  active 
material  is  assumed  to  expand  isotropically  upon  charging.  The 
battery  is  assumed  to  be  under  galvanostatic  operation.  Consid¬ 
ering  a  battery  in  which  the  cathode  and  anode  are  placed  face  to 
face,  lithium  ions  are  inserted  into  electrode  firstly  from  both  side 
faces.  If  the  effect  of  concentration  on  the  material  properties  is 
neglected  the  biaxial  modulus  and  partial  molar  volume  of  active 
plate,  £'  and  £?,  and  the  biaxial  modulus  of  current  collector,  E'c, 
then,  are  constant  and  independent  ofz-coordinate.  It  follows  from 
Eqs.  (14)-(16)  that 


The  closed-form  expressions  for  the  equal  biaxial  stress  in  the 
active  plates  and  current  collector  can  be  obtained, 


<Ti(Z,  t)  = 


inE'  £2Yi\ 
3  FD 


-  1  1_2  _ 

a  ~  Tz  + 


v(-ir 


7 r2  n2 

n= 1 


x  cos (httz)  exp(-n27T2t) 


(20) 


orc(z,  t)  = 


inE'c£2h i 
3  FD 


(1  -x)t 


(21) 


Substitution  of  Eq.  (9)  gives 

/hi  nhi  rhi 

E'dz  +  K  E'zdz=  -  E'£2c  dz 

hc-/i2  J -hc-/i2  J -hc-h2 

/hi  />hi  1  /*hi 

E'z  dz  +  K  /  E'z2  dz=  -  E'ficz  dz 

hc-h2  J -hc-h2  J -hc-h2 

Solving  for  £0  and  k  leads  to 

INC  -  BMC 

E°  ~  aTIp 

AMC  -  BNC 

K~  aTHp 


(12) 

(13) 


(14) 

(15) 


where  A,  B  and  I  are  biaxial  extensional  stiffness,  biaxial 
stretching-bending  coupling  stiffness  and  biaxial  bending  stiff¬ 
ness,  Nc  and  Mc  are  concentration  force  resultant  and  concentration 
moment  resultant,  given  by 


A  = 


rh, 

r*i\ 

r1 

/  E'  dz. 

B=  E'z  dz, 

i  = 

/ 

)  -hc -h2 

J  —hc~h  2 

j 

E'z 2  dz 


/hi  ph-L 

E'tic  dz,  Mc  =  1  /  i 

hc-h2  J -hc-h2 


E'Qcz  dz 


(16a-e) 


where  t  =  ,  z  =  /  =  (hc/h i  )/(2  +  (hc/hi  ){E'JE\ )) 

The  interesting  is  that  the  stress  in  the  current  collector  in 
the  electrode  does  not  vary  along  the  thickness  direction.  This  is 
because  the  plate  curvature  /c  is  equal  to  zero  in  the  large  and  sym¬ 
metric  layered  electrode  structure.  The  biaxial  stress  in  the  current 
collector,  expressed  in  Eq.  (9),  is  thus  irrelevant  with  the  coordinate 
z  in  the  thickness  direction. 

Fig.  2  shows  the  distribution  of  lithium  ions  concentration  along 
thickness.  As  lithium  ions  continuously  go  into  the  active  plate, 
the  lithiated  region  tends  to  expand,  while  the  unlithiated  region 
and  the  current  collector  tend  to  remain  undeformed.  Therefore, 
compressive  stress  is  generated  in  the  active  plates  due  to  the 
restriction  from  the  current  collector,  while  tensile  stress  is  caused 
in  the  current  collector  because  of  the  traction  of  active  plates.  Both 
stresses  keep  increasing  as  more  and  more  lithium  ions  migrate 
into  the  active  plates,  as  illustrated  in  Fig.  3.  In  the  calculation, 
the  current  collector  is  considered  as  copper,  and  the  active  plate 
is  assumed  to  be  graphite  plate.  Therefore,  the  modulus  ratio  of 
the  two  layers  F'JF\  is  chosen  to  be  10.  The  thickness  ratio  of  the 
two  layers  hc/hi  is  0.25,  determined  from  experimental  observation 
[12]. 

It  can  be  concluded  from  Figs.  2  and  3  that  the  stress  in  a 
large  active  plate  under  galvanostatic  charging  mainly  depends 
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Fig.  2.  Dimensionless  lithium  ions  concentration  c  (c  =  cFD/inh i)  along  the  plate 
thickness  under  galvanostatic  charging.  Here,  zjh\  =  0  represents  the  interface 
between  the  current  collector  and  active  plate,  z\h\  =  1  represents  the  side  face  of 
active  plate. 

on  Li  concentration.  The  region  with  higher  Li  concentration 
corresponds  to  a  larger  stress  level.  Therefore,  the  maximum  stress 
in  the  active  plate  is  found  to  be  at  the  side  face  where  lithium  ions 
firstly  migrate  into  the  active  plate. 

One  of  the  main  purposes  of  this  work  is  to  discuss  the  parame¬ 
ters  of  layered  structure  which  affect  the  diffusion  induced  stress. 
Here  are  two  important  parameters  from  mechanics  viewpoint,  i.e. 
the  modulus  ratio  of  current  collector  to  active  plate  as  well  as  the 
thickness  ratio  of  the  two  layers. 

Fig.  4  shows  the  variation  of  biaxial  stress  along  the  whole  elec¬ 
trode  thickness  with  respect  to  different  modulus  ratio  E'c/Ey  It 
clearly  shows  that  both  the  tensile  stress  in  current  collector  and 
the  compressive  stress  in  the  active  plate  reduce  as  the  ratio  E'JE\ 
declines,  but  more  remarkably  in  the  current  collector.  This  phe¬ 
nomenon  suggests  that  a  softer  and  more  flexible  current  collector 
would  help  to  reduce  the  restriction  between  the  two  layers,  and 
thereby  reducing  the  stress  in  the  whole  electrode.  Soft  material 
is  thus  preferred  in  the  material  selection  in  designing  of  current 


Fig.  3.  The  normalized  biaxial  stress  dr  (dr  =  cr-l/(E/1£2inh-i/3FD)  in  the  active  plate, 
and  dr  =  crc/(E^£2inh-i/3FD)  in  the  current  collector)  along  the  electrode  thickness 
under  galvanostatic  charging.  Here,  z/h\  <  0  represents  the  region  in  the  current 
collector,  z\h\  >0  represents  the  region  in  the  active  plate. 


Fig.  4.  Variation  of  the  normalized  biaxial  stress  o  along  thickness  with  different 
modulus  ratio  Fc/Fi  in  both  current  collector  and  active  plate. 

collector  on  the  condition  that  the  electrochemical  properties  are 
satisfied. 

Regarding  the  geometrical  configuration,  Figs.  5  and  6  show  the 
stress  in  the  current  collector  and  the  active  plate  with  respect  to 
different  thickness  ratio  hc/hi,  respectively.  As  shown  in  Fig.  5,  the 
layered  structure  electrodes  with  larger  thickness  ratio  hc/h \  result 
in  much  smaller  stress  in  the  current  collector  during  charging. 
Especially,  at  dimensionless  charging  time  t  =  0.8,  the  stress  in  cur¬ 
rent  collector  is  reduced  much  more  than  half.  Therefore,  a  larger 
thickness  ratio  is  beneficial  to  the  strength  of  current  collector. 

However,  thicker  current  collector  means  stronger  restriction  to 
the  active  plate.  Thinner  current  collector  could  elongate  more  dur¬ 
ing  charging,  to  some  extent  enabling  the  active  plate  more  relaxed. 
As  shown  in  Fig.  6,  thicker  current  collector  induces  larger  stress  in 
the  active  plate.  If  it  was  in  discharging,  the  thicker  current  collector 
would  induce  larger  tensile  stress  which  is  critical  to  mechanical 
failure  such  as  debonding  and  fracture. 

Therefore,  considering  the  strength  of  both  current  collector  and 
active  plate,  the  thickness  ratio  hc/fri  could  not  be  too  small  for 
safety  of  current  collector,  neither  could  it  be  too  large  in  order  to 
decrease  the  stress  in  the  active  plate  and  provide  enough  rooms 
for  the  active  materials.  From  mechanical  viewpoint,  the  current 


Fig.  5.  Variation  of  the  normalized  biaxial  stress  crc  ( dc  =  ocl(F\  C2inh\  / 3FD ))  in  the 
current  collector  with  respect  to  the  ratio  of  hc/h i  under  galvanostatic  charging. 
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Fig.  6.  Variation  of  the  normalized  biaxial  stress  b\  (di  =  0\  / (£'  £2inh\  /3 FD))  in  the 
active  plate  with  respect  to  the  ratio  of  hc//ii  under  galvanostatic  charging:  (a)  at 
the  interface  between  current  collector  and  active  plate,  z\h\  =  0;  (b)  at  the  side  face 
of  the  active  plate,  z\h\  =1. 


collector  should  be  as  thin  as  possible  on  the  precondition  that  the 
strength  is  satisfied. 


4.  Bilayer  electrode  plate 


Here  we  exam  the  bilayer  electrodes  where  a  single  active  plate 
is  bonded  to  current  collector.  By  substituting  h2  =  0  into  the  solu¬ 
tion  derived  in  Section  2  we  have  the  expressions  for  the  curvature 
and  the  stresses  in  the  active  layer  and  current  collector  as  follows, 
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Fig.  7.  Variation  of  dimensionless  plate  curvature  k  (k  =  KFD/inQ)  at  different 
charging  time  with  respect  to  the  thickness  ratio  hc/hi  in  bilayer  electrode. 
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In  bilayer  electrode,  the  most  important  feature  is  that  the  cur¬ 
vature  k  =f=  0.  As  result,  not  only  diffusion  induced  stress,  but  also 
bending  stress  will  evolve  during  charging. 

In  bilayer  electrode,  the  charging  process  is  the  same  as  that 
in  symmetric  electrode.  Fig.  7  shows  the  effects  of  thickness  ratio 
hc/hi  on  electrode  plate  curvature.  When  there  is  no  current  col¬ 
lector  (hc/hi  =0),  the  curvature  is  not  so  large.  The  deformation  of 
the  active  plate  is  purely  determined  by  Li  intercalation.  Once  a 
thin  current  collector  is  attached,  the  restriction  will  inhibit  the 
expansion  of  the  active  plate  in  the  region  adjacent  to  the  interface, 
and  thereby  increase  the  curvature.  Once  a  thicker  current  collector 
is  attached,  because  the  current  collector  usually  possesses  much 
higher  modulus,  its  large  flexural  rigidity  becomes  more  dominant, 
and  thus  reduces  the  overall  curvature  of  the  electrode.  This  is  the 
reason  why  the  curvature  k  is  found  to  be  firstly  promoted,  and 
then  reduced  with  respect  to  the  increase  of  hc/hi  in  Fig.  7. 

Fig.  8  shows  the  dimensionless  biaxial  stress  during  charging  in 
both  layers.  It  is  found  that  both  layers  subject  to  compressive  and 
tensile  stress  simultaneously.  Generally,  in  the  active  plate,  com¬ 
pressive  stress  is  induced  near  the  interface  because  the  lithiation 
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Fig.  8.  The  normalized  biaxial  stress  a  along  thickness  at  different  time  under  gal- 
vanostatic  charging.  Here,  z\h\  <  0  represents  the  region  in  the  current  collector, 
z\h\  >  0  represents  the  region  in  the  active  plate. 


induced  expansion  is  restricted  by  current  collector.  While  in  the 
region  near  the  side  face  the  stress  changes  from  compressive 
to  tensile  as  charging  goes  on.  This  is  the  result  of  competition 
between  diffusion  induced  stress  and  bending  stress.  At  beginning 
of  lithiation,  the  curvature  is  small.  The  overall  stress  is  mainly 
determined  by  diffusion  induced  stress.  As  more  and  more  lithium 
ions  migrate  into  the  active  plate,  the  curvature  increases  dramat¬ 
ically,  resulting  in  a  convex  curvature  at  the  side  face.  And  thus, 
the  stress  near  the  side  face  changes  to  tensile  because  bending 
becomes  more  important. 

Once  certain  amount  of  Li  has  been  intercalated,  the  peak  stress 
in  both  layers  is  found  at  the  interface.  It  suggests  that  considerably 
large  stress  near  the  interface  is  essential  to  make  the  deforma¬ 
tion  compatible.  Since  the  large  peak  stress  is  very  important  to 
structural  failure,  such  as  fracture,  here  is  a  key  question:  how  to 
reduce  the  interaction  between  layers,  and  thus  reduce  the  stress 
in  electrode  by  modifying  the  electrode  structure? 

Part  of  the  answer  can  be  found  in  Fig.  9,  which  shows  the  biax¬ 
ial  stress  in  the  current  collector  at  the  interface  upon  charging. 
The  stress  is  found  to  decrease  dramatically  with  increasing  of  cur¬ 
rent  collector  thickness  when  the  ratio  hc/h i  is  not  so  large.  This 


Fig.  9.  Variation  of  the  normalized  biaxial  stress  bc  at  the  interface  in  the  current 
collector  side  with  respect  to  the  ratio  of  hc/h-i  under  galvanostatic  charging. 


Fig.  10.  Variation  of  the  normalized  biaxial  stress  b\  in  the  active  plate  with  respect 
to  different  ratio  of  hc/hi  under  galvanostatic  charging:  (a)  at  the  interface  between 
current  collector  and  active  plate,  z\h\  =0;  (b)  at  the  side  face  of  the  active  plate, 
z/hi  =  1. 


decrease  may  be  attributed  to  the  increase  of  cross  section  area. 
However,  when  the  ratio  hc/hi  increases  continuously,  the  stress 
is  found  to  go  up  again.  This  may  be  attributed  to  the  decrease 
of  electrode  curvature  /c,  as  shown  in  Fig.  7.  When  the  ratio  hjh^ 
increases  from  around  0.1,  the  deformation  of  current  collector 
changes  from  bending-dominant  to  elongating-dominant,  and  thus 
slightly  increases  the  biaxial  stress.  On  the  whole,  increasing  the 
thickness  ratio  hc/hi  could  help  to  reduce  the  biaxial  stress  in  the 
current  collector. 

Fig.  1 0  shows  the  biaxial  stress  in  the  active  plate  with  respect  to 
different  thickness  ratio  hc/fri  upon  charging,  at  both  the  interface 
(z  =  0)  and  side  face  (z  =  1 ).  Unfortunately,  it  is  found  that  increas¬ 
ing  of  thickness  ratio  /ic/hi  would  increase  the  peaks  stresses,  no 
matter  compressive  or  tensile.  This  is  very  easy  to  understand  if 
we  draw  an  analogy  to  thermal  expansion.  It  is  the  restriction  of 
thermal  expansion,  not  the  expansion  itself,  induces  thermal  stress. 
Here,  increasing  of  thickness  ratio  represents  stronger  restriction. 
Therefore,  the  stresses  are  found  to  be  promoted  with  thicker  cur¬ 
rent  collector.  Again,  a  similar  conclusion  to  that  in  the  symmetric 
case  can  be  drawn  that  the  current  collector  should  be  as  thin  as 
possible. 
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Fig.  11.  Variation  of  the  normalized  biaxial  stress  o  along  electrode  thickness  with 
different  modulus  ratio  Ec/£i  in  the  bilayer  electrode. 


The  peak  stress  can  be  relaxed  by  modifying  the  material’s  elas¬ 
tic  property.  Fig.  1 1  shows  that  stress  in  the  current  collector  is 
greatly  reduced  by  decreasing  the  ratio  of  elastic  modulus  of  cur¬ 
rent  collector  to  active  plate,  i.e.  E'c/E'v  It  suggests  that  application 
of  appropriate  material  for  current  collector  could  help  to  reduce 
the  stress  in  both  layers.  As  discussed  above,  a  softer  and  more  flex¬ 
ible  current  collector  imposes  less  restriction  to  the  active  plate, 
and  therefore  reduces  the  stress  in  both  layers.  In  the  calculation 
shown  in  Fig.  11,  the  current  collector  is  thin  compared  with  the 
active  plate.  The  ratio  is  only  1 :4.  Therefore,  the  stress  reduction  is 
very  obvious  in  current  collector,  but  not  so  evident  in  active  plate. 


5.  Charging  conditions 

Another  purpose  of  this  work  is  to  address  the  effects  of  charging 
conditions.  In  literature,  most  works  implicitly  assume  that  lithium 
ions  migrate  into  all  active  particles  at  the  same  time  because  all 
particles  are  soaked  in  electrolyte.  However,  due  to  the  limited 
porosity  in  electrode  as  well  as  the  limited  diffusion  velocity  of 
lithium  ions  in  electrolyte,  it  is  very  difficult  for  all  active  particles 
to  react  with  lithium  ions  simultaneously.  In  experiment,  Li  concen¬ 
tration  gradient  in  active  plate  during  charging  has  been  observed 
[12].  It  is  necessary  to  investigate  the  effects  of  different  charging 
conditions. 

Influence  of  three  charging  conditions,  i.e.  uniform,  galvano- 
static  and  potentiostatic,  on  the  stress  in  active  plate  are  to  be 
discussed.  Due  to  the  different  charging  duration,  it  is  very  difficult 
to  compare  the  three  charging  process  with  respect  to  Newtonian 
time.  It  would  be  better  to  take  the  total  amount  of  intercalated 
lithium  ions  as  the  scale  of  charging  process.  Furthermore,  it  will 
be  revealed  from  following  derivation  that  the  expressions  for  the 
stresses  in  terms  of  the  total  amount  ions  are  much  simpler  than 
that  in  terms  of  Newtonian  time. 

For  simplicity,  the  symmetric  electrode  layered  structure  is  con¬ 
sidered.  The  influence  of  plate  curvature  is  thus  excluded.  Write  the 
dimensionless  total  amount  of  intercalated  lithium  ions  as 

ED  fhl 

c&dz  (25) 

Jo 

The  stress  during  uniform  charging  can  be  obtained  from  Eq.  (9): 
CTu’max  =  E\Qi2t/3FD  =  (26) 


Fig.  12.  The  maximum  diffusion  induced  stresses  d\  with  respect  to  dimensionless 
total  amount  of  intercalated  lithium  ions  in  three  charging  conditions:  (a)  overview 
in  which  the  potentiostatic  case  is  plotted  with  different  ratio  of  coFDIinh\ ;  (b) 
enlarged  view  to  show  the  intersection  of  stress  lines  corresponding  to  different 
charging  conditions. 


According  to  Fig.  3,  the  maximum  stress  in  non-uniform  charg¬ 
ing  is  found  at  the  side  face  where  z  =  1.  Combining  Eq.  (25)  with 
Eqs.  (20)  and  (23),  the  maximum  stress  can  be  obtained 
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for  galvanostatic  charging,  and 

ffl(l.a)  ,  n  c0FD 
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(27) 


(28) 


for  potentiostatic  charging.  Although  the  expressions  of  stress  in 
terms  of  Newtonian  time  is  very  complicated,  it  is  of  interesting 
that  they  are  very  simple  here.  The  stress  in  uniform  and  potentio¬ 
static  charging  is  exactly  linear  to  the  dimensionless  total  amount 
of  intercalated  lithium  ions  Q.  And  in  galvanostatic  case  with  suf¬ 
ficient  charging,  the  stress  is  also  linear  to  Q  with  expression 


=— XQ-1 
£'  Qinhy  /3FD  3 

since  the  progression  term  approaches  to  zero. 
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Fig.  12,  which  plots  the  stresses  against  the  dimensionless 
total  amount  of  intercalated  lithium  ions,  reveals  very  important 
features  of  the  maximum  diffusion  induced  stress  in  the  whole 
charging  history  with  respect  to  three  kinds  of  conditions.  Gen¬ 
erally,  the  maximum  stress  in  potentiostatic  charging  crp,  max  IS 
always  larger  than  that  in  uniform  case  (TU)max  during  whole  dura¬ 
tion  except  the  last  minute  when  the  electrode  is  saturated  and 

G  p,max  =  G  u.max- 

Comparing  uniform  and  galvanostatic  charging,  Fig.  12  shows 
that  the  stress  lines  of  the  two  stresses  become  parallel  after  suffi¬ 
cient  charging  with  a  dimensionless  difference  of  1  /3.  This  can  also 
be  verified  with  Eqs.  (26)  and  (29).  The  progression  term  in  Eq.  (27) 
only  makes  effects  in  the  early  stage  of  charging. 

Comparing  potentiostatic  and  galvanostatic  cases,  it  is  found 
that  there  is  a  cross  point  of  the  stress  lines  corresponding  to 
the  two  charging  conditions.  Combining  Eqs.  (27)  and  (29),  the 
cross  point  is  at  Q.  =  c0FD/inh i  -  1  /3.  It  means  that  (Tp.max  >  <Tg,max 
when  Q.  is  smaller  than  the  critical  value,  vice  verse,  <TP)max  <  <7g,max 
when  Q.  is  larger  than  the  critical  value.  This  suggests  an  optimized 
charging  operation  to  reduce  the  diffusion  induced  stress,  i.e.  first 
galvanostatic  charging  followed  by  potentiostatic. 

On  the  whole,  the  maximum  stress  in  potentiostatic  charging  is 
far  larger  than  the  other  two  cases  in  almost  whole  charging  his¬ 
tory.  The  maximum  stress  in  galvanostatic  charging  is  larger  than 
that  in  uniform  case  with  a  dimensionless  value  of  1  /3,  except  the 
beginning  stage.  Our  results  show  that  effects  of  Li  concentration 
gradient  on  stress  in  composite  level  cannot  be  neglected.  Calcula¬ 
tion  based  on  uniform  charging  would  result  in  smaller  peak  stress, 
and  in  turn,  gives  an  overestimation  on  electrode  strength. 

6.  Concluding  remarks 

In  this  paper,  layered  structures  applied  in  lithium  ion 
batteries  have  been  evaluated  by  formulating  the  diffusion 
induced  stress.  The  mechanical  role  of  current  collector  has  been 
addressed.  Finally,  different  charging  conditions  are  compared  and 
discussed. 

In  an  electrode  where  the  active  material  expands  upon  lithia- 
tion,  current  collector  acts  as  substrate  of  active  plate.  It  restricts  the 
expansion  of  active  plate,  and  also  deforms  with  the  plate.  In  a  large 
electrode  with  symmetric  layered  structure,  current  collector  and 
active  plates  are  elongated  without  bending.  Current  collector  is 
subjected  to  uniform  tensile  stress.  The  active  plates  are  subjected 
to  non-uniform  compressive  stress  which  heavily  depends  on  the 
distribution  of  Li  concentration.  The  peak  stress  in  the  active  plate 
is  found  at  the  side  faces  where  lithium  ions  firstly  migrate  into  the 
plate.  While  in  a  large  electrode  with  asymmetric  bilayer  structure, 


both  elongation  and  bending  are  found  during  lithiation.  Current 
collector  is  subjected  to  non-uniform  stress  caused  by  bending. 
Stress  in  the  active  plate  is  sum  of  bending  stress  and  diffusion 
induced  stress.  And  the  peak  stress  is  found  at  the  interface.  From 
viewpoint  of  mechanical  role,  thicker  current  collector  could  help 
to  reduce  the  stress  in  itself  while  increase  the  stress  in  active  plate. 
Therefore,  current  collector  should  be  as  thin  as  possible  based  on 
the  precondition  that  strength  is  satisfied.  It  is  also  found  the  cur¬ 
rent  collector  with  smaller  elastic  modulus  would  result  in  much 
lower  stress  in  electrode.  It  suggests  an  effective  way  to  release  the 
stress  in  electrode,  i.e.  choosing  softer  and  more  flexible  material 
for  current  collector.  Finally,  effects  of  three  charging  conditions,  i.e. 
uniform,  galvanostatic,  and  potentiostatic,  on  the  diffusion  induced 
stress  is  discussed.  It  is  found  the  maximum  stresses  for  three  cases 
are  linear  to  the  total  amount  of  intercalated  lithium  ions.  The  max¬ 
imum  stress  in  potentiostatic  charging  is  far  larger  than  the  other 
two  cases  in  almost  whole  charging  history.  The  maximum  stress 
in  galvanostatic  charging  is  larger  than  that  in  uniform  case  with 
a  dimensionless  value  of  1  / 3,  except  the  beginning  stage.  Based  on 
the  stresses,  an  optimized  charging  operation,  i.e.  first  galvanostatic 
followed  by  potentiostatic,  is  suggested. 
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